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Preface 


In  late  1976,  a  study  to  produce  a  wave  climate  for  U.  S.  coastal 
waters  was  initiated  at  the  U.  S.  Army  Engineer  Waterways  Experiment 
Station  (WES).  The  Wave  Information  Study  (WIS)  was  authorized  by  the 
Office,  Chief  of  Engineers,  U.  S.  Army,  as  a  part  of  the  Coastal  Field 
Data  Collection  Program  which  is  managed  by  the  U.  S.  Amy  Coastal  En¬ 
gineering  Research  Center.  The  U.  S.  Army  Engineer  Division,  South 
Atlantic,  and  the  U.  S.  Amy  Engineer  Division,  New  England,  also 
authorized  funds  during  the  initial  year  of  this  study  (FY  1978)  to 
expedite  execution  of  the  Atlantic  coast  portion  of  this  program. 

This  report,  the  eighth  in  a  series,  presents  the  theoretical  de¬ 
velopment  of  the  Phase  III  portion  of  WIS  which  is  employed  to  calculate 
wave  information  in  shallow  water.  The  study  was  conducted  in  the  Hy¬ 
draulics  Laboratory  under  the  direction  of  Mr.  H.  B.  Simmons,  Chief  of 
the  Hydraulics  Laboratory,  Dr.  R.  W.  Whalin,  fomer  Chief  of  the  Wave 
Dynamics  Division,  and  Mr.  C.  E.  Chatham,  Jr.,  Acting  Chief  of  the  Wave 
Dynamics  Division.  Mr.  J.  H.  Lockhart,  Jr.,  is  the  Office,  Chief  of  En¬ 
gineers,  Technical  Monitor  for  the  Coastal  Field  Data  Collection  Program 
This  report  was  prepared  by  Dr.  R.  E.  Jensen  with  assistance  from 
Mr.  W.  D.  Corson  and  Ms.  J.  H.  Messing.  Initial  work  in  the  theoretical 
treatment  of  the  Phase  III  methodology  was  performed  with  the  help  of 
Dr.  D.  T.  Resio. 

Commanders  and  Directors  of  WES  during  the  conduct  of  the  study 
and  the  preparation  and  publication  of  this  report  were  COL  John  L. 
Cannon,  CE,  COL  Nelson  P.  Conover,  CE,  and  COL  Til  ford  C.  Creel,  CE. 
Technical  Director  was  Mr.  F.  R.  Brown. 
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METHODOLOGY  FOR  THE  CALCULATION  OF  A 
SHALLOW-WATER  WAVE  CLIMATE 


Introduction 


1.  The  primary  purpose  of  the  Wave  Information  Study  (WIS)  is  to 
produce  useful,  reliable  wave  data  for  planning  and  design  in  coastal 
waters.  During  the  last  6  years,  a  methodology  has  evolved  to  simulate 
the  offshore  wave  climate  via  numerical  procedures.  Extensive  efforts 
have  been  directed  into  problems  of  wind  and  wave  prediction  over  large 
ocean  areas  from  historical  records.  Reports  1-7  in  this  series  docu¬ 
ment  various  portions  of  this  work.  From  comparisons  of  hindcast  and 
measured  data  sets,  it  would  appear  that  the  climatological  character¬ 
istics  of  the  simulated  data  in  the  offshore  area  are  quite  similar  to 
those  of  the  measured  data  (Corson  and  Resio  1981).  The  task  of  accu¬ 
rately  transforming  these  offshore  data  into  the  nearshore  region  now 
assumes  a  critical  role  in  the  final  stage  of  obtaining  the  overall  WIS 
goal. 

2.  Previous  approaches  to  nearshore  wave  transformations  have 
concentrated  on  the  calculation  of  a  limited  number  of  wave  conditions 
over  very  small  spatial  areas.  The  reason  for  this  is  that,  typically, 
the  results  are  needed  at  a  particular  site  for  some  specific  design  or 
planning  consideration.  While  this  approach  is  convenient  when  trans¬ 
forming  a  small  number  of  given  offshore  wave  conditions  into  nearshore 
wave  conditions,  it  becomes  very  cumbersome  when  attempting  to  define  a 
nearshore  wave  climate  over  20  years  time  along  thousands  of  miles  of 
coast.  Instead,  the  WIS  approach  has  been  to  analyze  successive  nested 
scales  of  natural  organization  in  time  and  space  along  a  coast.  As 
shown  in  Figure  1,  this  leads  to  a  three-phase  organization  of  WIS,  with 
various  processes  becoming  dominant  at  different  distances  from  a  coast. 

3.  It  is  impossible  to  foresee  every  site  at  which  detailed  wave 
information  will  be  required  for  specific  projects  in  the  future.  Also, 
it  is  not  cost-effective  to  perform  detailed  calculations  on  an  ultra- 
fine  grid  for  all  coastal  areas,  since  most  of  this  information  would 
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Figure  1.  Schematic  representation  of  the  three 
phases  of  the  WIS 


never  be  needed.  Consequently,  the  Phase  III  data  must  be  computed  in 
such  a  manner  as  to  be  representative  of  some  length  of  coast  which  is 
sufficient  to  be  useful  for  most  projects,  while  sufficiently  large  to 
filter  out  many  spurious  local  effects.  As  shown  in  Plates  1-6,  10-mile 
coastal  segments  were  selected  to  meet  this  requirement.  Small-scale 
variations  in  wave  climate  within  a  10-raile  segment  can  be  treated  as 
perturbations  on  this  Phase  III  scale  if  additional  spatial  detail  is 
required. 

A.  Several  advantages  are  associated  with  the  selection  of  this 
10-mile  scale  in  Phase  III.  The  most  notable  of  these  is  probably  the 
dramatic  decrease  in  cost  of  computation.  This  is  particularly  important 


since  the  processes  involved  in  nearshore  wave  transformation  are  not 
completely  understood  at  the  present  time  (Vincent  1982).  It  is  ex¬ 
pected  that  results  in  the  nearshore  zone  may  have  to  be  recomputed  at 
some  later  time  as  an  understanding  of,  and  ability  to  model,  wave 
transformation  processes  increases.  Thus  it  would  not  seem  reasonable 
to  develop  extremely  costly,  time-consuming  methods  for  the  present 
Phase  III  calculations. 

5.  Although  the  techniques  used  must  be  general  in  terms  of  dis¬ 
tance  along  the  coast,  recent  field  evidence  has  suggested  that  signifi¬ 
cant  variations  in  wave  energy  occur  systematically  in  the  onshore- 
offshore  direction.  These  effects  will  not  average  out  over  longshore 
distance;  consequently,  all  Phase  III  results  must  be  referenced  to  a 
particular  depth.  In  very  shallow  water,  the  effects  of  depth  become  so 
dominant  that  variations  in  mean  water  levels  (tides,  surges,  setup, 
etc.)  are  the  controlling  factors  in  the  external  portion  of  the  wave- 
height  distribution.  Three  situations  in  wave  climate  estimations  can 
be  related  to  the  relative  importance  of  water-level  variations.  In  the 
innermost  and  outermost  regions,  water-level  variations  and  wave-height 
variations  can  be  treated  independently.  However,  in  the  intermediate 
region,  irfilch  typically  is  the  area  of  primary  interest,  the  benavior  of 
both  must  be  considered  together  in  the  formulation  of  a  wave  climate. 
For  this  first  Phase  III  report,  a  10-m  depth  has  been  chosen  as  repre¬ 
sentative  of  nearshore  conditions  while  not  being  too  dependent  on  most 
water-level  variations. 

6.  It  should  also  be  noted  that  the  present  report,  like  the  pre¬ 
vious  Phase  I  and  Phase  II  reports,  is  limited  to  extratropical  condi¬ 
tions.  That  portion  of  the  wave  climate  contributed  by  tropical  storms 
will  be  added  to  this  data  set  in  a  later  report. 

Development  of  Phase  III  Methodology 


Brief  perspective  of 
previous  methodologies 

7.  Many  different  techniques  are  presently  used  to  transform 


deepwater  wave  information  into  shallow  water.  In  general,  these 
techniques  can  be  separated  into  two  broad  classes,  those  which  trans¬ 
form  individual  wave  conditions  and  those  which  transform  probability 
distributions  of  wave  parameters  from  one  space  (or  water  depth)  to 
another. 

8.  The  first  class  of  techniques  generally  transforms  a  specified 
wave  condition  at  a  boundary  (defined  by  a  set  of  wave  parameters), 
through  some  intermediate  water  depth  region,  into  shallow  water.  Typi¬ 
cally,  this  method  is  based  on  linearized  formulations  of  monochromatic 
wave  theories  where  the  dependent  wave  parameters  are  influenced  by 
shallow-water  effects  (wave  shoaling,  refraction,  diffraction,  bottom 
friction,  breaking,  etc.)  (Skovgaard,  Jonssen,  and  Bertelsen  1975). 

More  recently,  computer  codes  based  on  approximations  to  dispersive  long 
waves  have  been  used  to  calculate  shallow-water  wave  transformations 
which  include  additional  nonlinear  effects.  However,  these  codes  tend 
to  be  quite  expensive  to  run;  and  in  most  open-coast  areas,  the  results 
obtained  are  quite  similar  to  linear  theory.  A  third  approach,  used  in 
some  studies,  involves  the  transformation  of  wave  spectra  propagating 
into  shallow  water.  This  method  usually  requires  extensive  amounts  of 
computer  time  since  each  frequency-direction  component  of  the  spectrum 
must  be  propagated  and  have  its  energy  source/sink  integration  performed 
individually  (Hsiao  1978). 

9.  In  the  second  class  of  techniques,  a  Jacobian  matrix  operator 
is  used  to  transform  distribution  of  wave  heights,  periods,  and  direc¬ 
tions  in  deep  water  into  corresponding  shallow-water  distributions  as 
shown  in  Appendix  A. 

10.  A  major  problem  in  this  approach  lies  in  the  evaluation  of 
the  Jacobian.  If  an  analytical  method  is  applied  for  this  purpose,  only 
certain  classes  of  functions  can  be  admitted  into  the  interrelationships 
between  different  wave  parameters.  This  restriction  could  be  overcome 
through  a  purely  numerical  approach;  however,  the  method  would  still 
lose  all  information  on  specific  events  (i.e.  individual  storms)  and 
would  not  provide  a  very  accurate  description  of  the  external  portion 
of  the  nearshore  wave  population. 


11.  Since  neither  of  the  traditional  approaches  to  the  estimation 
of  nearshore  wave  characteristics  seemed  appropriate,  a  third  approach 
to  the  problem  of  nearshore  wave  transformations  was  developed  specifi¬ 
cally  to  meet  the  needs  of  the  Phase  Ill  portion  of  WIS.  It  should  be 
recognized  at  the  outset  of  this  discusssion  that  the  primary  objective 
of  this  new  approach  is  to  maximize  calculation  efficiency  rather  than 
to  introduce  a  new  wave  transformation  technique.  This  increase  in  ef¬ 
ficiency  is  accomplished  by  recognizing  that  over  sufficiently  small  in¬ 
tervals,  wave  parameters  may  be  considered  approximately  equal.  In 
other  words,  a  1.01-m  wave  height  is  not  essentially  different  from  a 
1.00-m  height  in  most  applications,  an  8.42-sec  wave  period  is  approxi¬ 
mately  equivalent  to  an  8.4-sec  wave  period,  and  so  on. 

12.  The  first  problem  which  must  be  resolved  in  this  approach  is 
how  many  wave  parameters  are  required  to  represent  deepwater  wave  condi¬ 
tions  such  that  all  basic  information  affecting  the  nearshore  transfor¬ 
mation  process  is  retained.  If  this  number  is  sufficiently  small,  and 
if  the  number  of  increments  needed  to  cover  the  range  of  parameter 
values  expected  is  reasonable,  then  a  discretized  approach  is  quite  ad¬ 
vantageous  to  performing  individual  calculations.  For  example,  if  the 
computations  were  made  on  an  individual  basis  at  each  Phase  III  site, 
almost  10  million  transformations  would  be  required  to  span  the  20-year 
period  of  the  hindcast  at  3-hr  intervals.  However,  if  only  three  param¬ 
eters,  each  spanned  by  20  intervals,  were  required  to  obtain  adequate 
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resolution,  then  only  (20)  or  8,000  transformation  calculations  would 
be  required.  More  generally,  an  n-parameter  transformation  would  re¬ 
quire  N*  calculations, 

n 

N  =  /7  m,  (I 

i=l 

trtiere  m^  is  the  number  of  intervals  used  to  resolve  the  i**1 


For  convenience,  symbols  are  listed  and  defined  in  the  Notation 
(Appendix  D) . 


parameter.  Since  these  transformations  have  been  generalized  for  paral¬ 
lel  bottom  contours,  these  calculations  would  suffice  for  all  166  shore¬ 
line  segments  for  a  given  depth. 

13.  It  is  important  to  note  here  that  other  than  the  parallel 
bottom  assumption,  no  assumptions  have  been  made  regarding  the  nature  of 
the  transformation  process.  Consequently,  any  number  of  transformation 
mechanisms  can  be  considered  within  the  framework  of  this  approach. 

This  means  that  it  is  easy  to  add  to  or  delete  mechanisms  from  the  com¬ 


puter  analog.  This  algorithm  is  used  to  establish  a  matrix  of  N  trans 

.  th 


formed  values  for  each  of  the  i  parameters.  When  the  individual 
transformations  are  performed,  the  incremental  deepwater  parameter 
values  can  then  be  viewed  as  indices  to  a  matrix  of  transformed  values; 
thus  we  have  for  any  individual  n-parameter  transformation  as  a  matrix 
formulation  given  by 


Pi  =  V1!’  V 


’  V 


p2  =  T2(ir  V  V 


(2) 


P'  =  T  (I 


n 


1’  2’ 


where  P!  are  the  nearshore  wave  parameter  values,  T^  are  the  trans¬ 


formation  matrices,  and  I  are  the  incremental  values  of  the  deepwater 


parameters.  The  only  calculations  required  are  the  conversions  of  the 


continuous  values  of  the  deepwater  parameters  P^  into  the  referencing 


indices 


1^  ,  which  usually  require  only  a  single  multiplication. 


Representation  of  Offshore  Wave  Conditions 
for  Wave  Transformation 


14.  In  the  Phase  II  wave  model,  20  discretized  frequency  incre¬ 
ments  and  16  discretized  angle  increments  were  used  to  approximate  con¬ 
tinuous  two-dimensional  wave  spectra  (Resio  1981).  If  all  320  elements 
in  the  16-X-20  array  were  treated  as  independent,  the  value  of  N  would 
become  ridiculously  large  in  an  analysis  of  the  type  just  described.  On 
the  other  hand,  if  the  spectra  were  reduced  to  simple  monochromatic  wave 


trains,  significant  information  important  to  nearshore  wave  transforma¬ 
tions  would  be  lost.  An  approach  that  attempts  to  retain  the  organized 
structure  of  the  spectra  and  to  include  these  effects  into  nearshore 
wave  transformations  would  be  to  use  a  parametric  representation  of  the 
offshore  wave  climate.  The  solutions,  although  they  can  be  computed  as 
two-dimensional  spectra,  are  effectively  limited  to  the  parameter  space 
itself.  Hasselmann  et  al.  (1973  and  1976)  give  good  descriptions  of  the 
use  of  parametric  models  for  wave  generation  and  propagation.  The 
principle  in  this  case  is  quite  similar. 

15.  Most  parameterizations  of  wave  spectra  to  date  have  treated 
the  frequency  variation  and  angular  variation  of  energy  density  sepa¬ 
rately.  Thus  the  two-dimensional  spectrum  can  be  written  as 

E2(f,6)  =  E1(f)<K6) 
or 

E2(f,0)  =  E1(f)<t>(0,f) 

where  E2  is  the  two-dimensional  spectrum  (in  frequency  f  and  direc¬ 
tion  0  ),  <]>  is  a  normalized  angular  distribution  function  and  _s 

the  one-dimensional  spectrum  defined  as 

2  71 

Ej (f )  =  J  E2(f,0)d0 
0 

The  number  of  parameters  used  to  represent  a  spectrum  varies  as  a  func¬ 
tion  of  the  detail  one  wishes  to  retain.  The  classic  Pierson-Moskowitz 
spectrum  (Pierson  and  Moskowitz  1964)  can  be  written  in  the  form 


where  g  is  acceleration  due  to  gravity,  f  is  frequency,  and  fm  is 
the  frequency  of  the  spectral  peak.  Essentially,  there  is  only  one  free 


parameter,  f  ,  in  this  representation.  More  recently,  the  JONSWAP 
spectrum  containing  three  free  parameters  was  introduced  by  Hasselmann 
et  al.  (1973) 

^ = & exp  ['HkTY 

where  a  ,  f  ,  and  6  are  all  free  parameters  and  v  is  a  given  func 

tion  of  (f/f  ).  The  JONSWAP  spectrum  is  somewhat  cumbersome  to  treat 
m 

analytically  and  for  all  but  the  last  stage  of  wave  growth  v  is  taken 
to  be  a  constant  equal  to  3.3.  In  many  recent  field  experiments  of 
active  wave  growth,  the  JONSWAP  spectrum  has  been  shown  to  provide  a 
reasonably  good  approximation  to  the  measured  spectrum;  however,  an  al¬ 
most  equivalent  but  much  simpler  form  for  these  spectra  was  postulated 
by  Kitaigordskii  (1962)  on  the  basis  of  similarity  theory 


This  form  has  only  two  free  parameters  a  and  f^  ;  consequently,  a 
simple  relationship  between  total  wave  energy  and  these  parameters  exists 
i.e. 

00 

E Q  =/  Ej(f)df  =  AorgV4  (6) 

where  \  is  a  constant. 

16.  Equation  5  should  provide  a  good  description  of  the  one¬ 
dimensional  spectrum  under  active  wave  generation  conditions;  therefore 
it  provides  a  description  of  what  is  typically  referred  to  as  the  local 
sea.  The  spectral  shape  could  be  made  increasingly  complex  by  introduc¬ 
ing  additional  spectral  parameters.  Since  the  portion  of  the  spectrum 
occupied  by  the  local  sea  has  been  defined,  the  major  contributions  re¬ 
maining  in  the  spectrum  must  come  from  swell  advected  into  coastal  areas. 


Since  the  parameters  a  and  f  are  contained  within  all  Phase  II  wave 

m 

records,  the  distinction  between  sea  and  swell  here  is  not  at  all 
arbitrary. 

17.  This  separation  of  the  wave  records  into  sea  and  swell  per¬ 
mits  distinctly  different  approaches  to  the  treatment  of  the  frequency 
and  directional  characteristics  of  these  wave  trains.  Swell  that  has 
propagated  beyond  its  area  of  generation  usually  contains  almost  all  of 
its  energy  in  a  narrow  frequency  and  direction  band.  Thus  swell  can  be 
adequately  described  as  a  unidirectional,  monochromatic  wave  train 
impinging  on  a  coast.  On  the  other  hand,  the  local  sea  tends  to  be 
quite  broad-banded  in  direction  as  well  as  frequency.  Therefore  a  di¬ 
rectional  distribution  function  must  be  included  for  the  local  sea  in 
order  to  account  for  this.  Resio  (1978)  has  shown  that  longshore  energy 
fluxes  (and  associated  longshore  sediment  transport  rates)  can  be  very 
significantly  affected  by  this  directional  distribution  of  energy  in  a 
spectrum;  thus  it  is  important  to  include  such  effects  in  the  nearshore 

transformation.  A  reasonable  approximation  to  average  directional  dis- 

4 

tributions  for  observed  spectra  is  given  by  a  cos  8  distribution  of 
energy  for  deepwater  input  conditions  to  the  Phase  III  parametric  spec¬ 
tral  form 

E2(f,0)  =  E^f)  cos4  (6  -  6)  (7) 


where  E^(f)  is  given  by  Equation  5  and  0  is  the  central  angle  of  the 
spectrum. 

18.  As  additional  information  accumulates  relative  to  a  better 
definition  of  directional  spectral  characteristics  (Vincent  1982),  im¬ 
proved  directional  distribution  functions  can  be  inserted  into  input 
conditions  to  the  Phase  III  data  if  desired. 


19.  At  this  point,  each  of  the  two  kinds  of  wave  trains  seen 
along  a  coast  have  been  shown  to  be  representable  by  three  parameters: 
wave  height  (H) ,  wave  period  (T) ,  and  wave  direction  (8)  for  the  unidi 


rectional  monochromatic  swell  and  a  ,  f  ,  and  8  for  the  local  sea. 

’  m 


For  consistency,  if  we  define  a  period  of  the  spectral  peak  as 


T  =  f 
P 


-1 

m 


and  use  Equation  6  to  obtain  wave  height  as  a  function  of  a  and  T 

p 

then  both  wave  trains  can  be  described  by  comparable  parameters  H  ,  T 
and  0  .  It  should  be  recognized  that  this  by  no  means  restricts  the 
local  sea  to  unidirectional  monochromatic  waves,  but  rather  only  to  a 
fixed  deepwater,  spectral  shape. 

20.  Due  to  the  nature  of  the  definition  of  local  sea,  only  one 
such  wave  train  can  exist  in  one  place  at  one  time;  however,  this  is  not 
the  case  for  swell,  since  it  can  propagate  from  various  source  regions 
into  a  particular  coastal  area.  An  analysis  of  wave  generation  areas  in 
the  Atlantic  Ocean  indicates  that  one  swell  wave  train  should  suffice  to 
provide  an  adequate  description  of  the  significant  aspects  of  nearshore 
wave  climate  generated  by  storms  which  have:  (a)  passed  a  site  and  are 
moving  away  or  (b)  passed  along  the  coast  some  distance  offshore  (Fig- 


21.  Figure  3  summarizes  the  previous  discussion  regarding  the 
representation  of  wave  propagation  into  the  coast.  As  seen  here,  two 
independent  wave  trains  are  considered  superimposed  at  the  coast:  a  lo¬ 
cal  sea  described  by  three  free  parameters  defining  a  two-dimensional 
spectrum,  and  a  swell  wave  train  from  passing  storms  described  by  three 
free  parameters  defining  unidirectional,  monochromatic  waves. 


SWELL  TRAIN  1 

INITIAL  CONDITIONS: 
-MONOCHROMATIC 
-UNIDIRECTIONAL 


WIND -GENERATED  SEA 

INITIAL  CONDITIONS: 

— 2-D  DISCRETE  SPECTRUM 
(KITAIGOR  DSKII-TYPE) 
-ASSUMED  SYMMETRIC  ABOUT 
CENTRAL  ANGLE 


Figure  3.  Generalized  deepwater  input  conditions  to 
the  Phase  III  transformation  technique 


Processes  Involved  in  Nearshore  Wave  Transformations 


22.  All  of  the  processes  that  act  on  wave  conditions  during  their 
transformation  from  deep  to  shallow  water  may  be  classified  as  to 
whether  they  are  conservative  or  nonconservative.  In  conservative  pro¬ 
cesses,  any  variations  in  spectral  energies  are  due  to  redistribution  of 
energy  density  rather  than  a  loss  (or  gain)  of  energy  out  of  (or  into) 
the  system.  Except  for  continued  wind  input,  all  of  the  nonconservative 
processes  for  shoaling  waves  represent  net  sinks  of  energy.  Thus,  in 
almost  all  cases,  there  is  a  net  loss  of  wave  energy  in  waves  propagat¬ 
ing  into  shallow  water. 


Conservative  processes 

23.  Traditionally,  the  primary  conservative  processes  considered 
in  nearshore  wave  transformation  have  been  refraction  and  shoaling. 
Models  have  been  developed  to  treat  these  processes  similar  to  the 
principles  of  geometric  optics.  Along  a  wave  ray,  the  fundamental  con¬ 
servation  law  is  given  by 


F^(k)  =  constant  (8) 

where  is  the  two-dimensional  wave  spectrum  in  wave  number  (k)  space 

Transforming  Equation  8  into  frequency-direction  space  yields 

E2(f,8)ccg  =  jE2(f,6)ccJ  (9) 

where  c  and  c^  are  the  phase  velocity  and  group  velocity  of  the  wave 
component,  respectively,  and  the  subscript  "o"  denotes  deepwater  condi 
tions.  Equation  9  states  that  in  this  space,  the  value  of  E2(f,0)cc^ 
is  always  equal  to  its  initial  (deep  water)  value.  Hence,  as  c  and  c 
change,  E2(f,0)  must  also  change. 

24.  Using  a  straight  parallel  bottom  assumption,  a  simple  spec¬ 
tral  model  for  refraction  and  shoaling  can  be  developed.  Since  we  are 
treating  a  relatively  broad  spectrum  in  the  case  of  the  local  sea,  lin¬ 
ear  wave  theory  should  be  quite  adequate  to  model  the  various  individual 
components.  Hence,  c  and  c^  are  functions  only  of  frequency  and 
depth.  This  gives  for  the  transformation  of  a  given  component 

E2(f,0)  =  p(f,h)|E2(f,0)j  (10) 

where  p(f,h)  is  defined  as: 

cc  \ 

_!)o 
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P(f,h)  = 


25.  It  can  be  seen  from  Equation  10  that  the  energy  density  trans 
formations  are  independent  of  approach  angle.  This  may  seem  strange  at 


first  since  the  refraction  process  is  actually  quite  sensitive  to  ini¬ 
tial  wave  angle.  This  paradox  is  resolved  when  a  representation  of  the 
total  energy  is  examined. 


®  In 

Et  =  J  J  E2(f,0)d0df  (11) 

0  0 


which  can  be  written  in  terms  of  the  initial  spectrum  as 


'll " 


2(f,0)]op(f,h)  deodfo 
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where  0^  is  used  to  denote  the  deepwater  wave  angle.  This  follows 
from  Equation  11  since  the  frequency  along  a  ray  remains  constant 
whereas  the  angle  does  not.  The  actual  working  limits  of  the  integra¬ 
tion  about  the  directional  distribution  of  E„(f,0  )  are  ±n/2  about 

2  1  o 

the  principal  wave  angle  propagation  normal  to  a  given  shoreline  segment 
All  energy  "outside"  of  this  region  propagates  seaward  and  thus  has  no 
effect  on  the  nearshore-wave  climate. 

26.  From  Snell's  law,  the  shallow-water  angle  can  be  calculated 
as  a  function  of  frequency,  depth,  and  the  initial  angle  (Appendix  B) . 


c  sin  0 

0(f,d,0  )  =  sin  — - - 

o  c 


(13) 


It  can  be  seen  from  Equations  12  and  13  that  the  effect  of  refraction  is 
to  compress  the  energy  density  into  a  successively  narrower  angle  band 
as  waves  propagate  toward  a  coast.  Even  though  the  energy  density  will 
remain  high,  the  phase  volume  can  significantly  decrease  for  certain  ap¬ 
proach  angles  and,  consequently,  so  can  the  total  energy. 

27.  A  numerical  treatment  of  Equation  12  is  very  simple  since 

00  I 

p(f,h)  and  are  independent  of  energy  density.  They  need  be 

do 

computed  only  once  for  a  given  set  of  frequency-direction  increments 
and  the  total  wave  energy  can  be  quickly  calculated  by  an  inner  product 
formulation 


where  i  and  j  represent  frequency  and  angle  increments,  respectively, 


and  where  A. .  represents  a  matrix  of  precomputed  values  of  p(f ,h) 


ij 
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multiplied  by  r^r— 

Oo 


28.  A  second  kind  of  conservative  process  consists  of  wave-wave 
interactions.  These  interactions  are  constrained  such  that: 


If 


3F2(k) 


at 


dk  =  0 


(15) 
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3n(k) 


3t 


dk  =  0 


(16) 
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SEj (f ) 


at 


df  =  0 


(17) 


where  n(k)  is  action  density  defined  in  terms  of  the  wave  number 
modulus  k  .  Equations  15,  16,  and  17  represent  conservation  of  total 
momentum,  total  action,  and  total  energy,  respectively.  These  equations 
do  not  constrain  these  quantities  to  remain  constant  for  a  particular 
frequency  direction  wave  component,  however,  since  energy,  momentum,  and 
action  can  be  interchanged  among  different  components.  Only  the  net  in¬ 
tegrals  of  the  exchanges  must  be  zero. 

Nonconservative  processes 

29.  Most  theories  of  energy  loss  in  shoaling  waves  have  centered 
on  wave-bottom  interactions  of  a  dissipative  nature.  Bottom  friction, 
percolation,  and  viscoelastic  effects  have  all  been  shown  theoretically 
to  be  capable  of  removing  significant  amounts  of  energy  from  waves  in 
shallow  water  (Hsiao  1978).  Unfortunately,  the  empirical  coefficients 
that  are  necessary  when  using  these  theories  are  quite  dependent  on 
sedimentary  characteristics  of  the  bottom  over  which  the  waves  pass. 

This  sensitivity  would  seem  to  suggest  that  significant  variation  of 
energy  loss  in  shoaling  waves  might  occur  over  relatively  small  time  and 


space  scales.  Furthermore,  since  these  characteristics  can  change 
through  time,  no  single  transformation  might  be  applicable  at  one  site 
for  all  time. 

30.  In  two  recent  studies  (Goda  1974  and  Battjes  1972),  irregular 
wave  breaking  has  been  used  to  attempt  to  derive  a  rate  of  energy  loss 
for  waves  near  a  coast.  In  both  approaches,  the  probability  of  wave 
height  exceeding  some  critical  value  is  assumed  to  be  governed  by  a 
modified  Rayleigh  distribution;  however,  what  happens  to  the  "broken" 
energy  is  treated  differently  in  the  two  approaches.  The  observation 
of  increases  in  whitecapping  in  shoaling  waves  and  extensive  surf  zones 
with  almost  continuous  wave  breaking  extending  well  offshore  during 
storms  are  two  supporting,  albeit  qualitative,  piece*-  of  evidence  for 
this  kind  of  wave  breaking.  The  problem  of  estimating  the  energy  loss 
described  by  Goda  and  Battjes  is  somewhat  analogous  to  the  problem  of 
estimating  energy  losses  in  turbulence  by  following  an  individual  eddy. 
To  begin  with,  like  the  eddies,  waves  are  a  superposition  of  many  inter¬ 
acting  components;  hence,  if  wave  breaking  occurs,  there  is  considerable 
ambiguity  as  to  which  components  are  breaking.  This  closure  problem 
makes  it  difficult  to  use  the  energy  loss  rates  based  on  this  approach 
with  much  confidence. 

31.  Recent  studies  by  Heterich  and  Hasselmann  (1980)  and  Shemdin 
et  al.  (1980)  have  shown  that  wave-wave  interactions  can  create  a  strong 
flux  of  energy  into  high  frequencies  in  shoaling  waves.  Although  the 
wave-wave  interactions  are  themselves  conservative,  this  flux  is  into  a 
region  where  the  waves  are  unstable  and  wave  breaking  (or  high  viscous 
losses)  can  occur.  The  advantage  of  this  approach  to  wave  breaking  is 
that  the  quantity  predicted  is  the  flux  of  energy  into  this  unstable 
region  and  hence  no  information  is  required  regarding  the  nature  of  the 
partitioning  of  "broken"  and  "unbroken"  energy  in  individual  waves.  Nor 
is  it  necessary  to  reformulate  a  revised  probability  law  of  individual 
wave  heights  in  shallow  water. 

32.  The  dominant  transformation  process  in  shallow  water  appears 
consistently  to  produce  a  similarity  form  for  a  spectrum  (Kitaigordskii , 
Krasitskii,  and  Zaslavaskii  1975).  In  deep  water,  as  seen  in 


Equations  3,  4,  and  5,  the  wave  spectrum  follows  an  f  spectral  shape 
in  the  region  greater  than  f  .  In  finite  water  depths,  the  general 
form  is  given  by 


a*V5 

Et(f)  =  -8-  -t  ♦(«.) 

1  <2n)4  ^ 


where  ♦(iu^)  is  a  nondimensional  function  varying  from  1.0  in  deep 
water  to  0.0  when  h  =  0.0  ,  as  shown  by  Figure  4.  When  uj^  is  less 
than  1.0,  <Kiuh)  can  be  approximated  by: 

*(V  -  |  (19) 

where 


F{  (jj  <h/g),/4  ) 

Figure  4.  The  functional  relationships  between  ♦(ui,  )  (solid  line) 

2  1/2 

and  ul/2  versus  (h/g)  '  (from  Kitaigordskii ,  Krasitskii,  and 

Zaslavaskii  1975) 

Replacing  Equation  20  into  Equation  18  will  modify  the  spectral  form  to 
an  f  shape  or: 


V(f  } 
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In  shallow  water,  the  limiting  spectral  form  is  related  to  the  water 

-3 

depth  h  and  is  proportional  to  an  f  spectral  shape.  Studies  by 
Goda  (1974),  Thornton  (1977),  Ou  (1980),  Iwata  (1980),  and  Vincent 
(1981)  support  the  form  given  by  Equation  21.  Employing  Equation  21  as 
the  shallow-water  spectral  shape  assumes  that  a  remains  constant  from 
the  wave  conditions  found  in  deep  water.  What  has  been  found  in  proto¬ 
type  data  is  that  a  varies  with  the  "significant  steepness"  of  the 
waves  (based  on  Hg)  or  the  dimensionless  peak  wave  number.  Since  bottom 
friction  effects  are  ignored  in  the  transformation  effects,  the  assump¬ 
tion  of  a  constant  a  is  not  a  bad  approximation.* 

33.  Since  the  geometry  of  the  wave-wave  interactions  appears  to 
force  a  constant  spectral  shape  in  wave  number  space,  the  ratio  of  mo¬ 
mentum  fluxes  in  various  portions  of  the  spectrum  retains  about  the  same 
relative  values  in  both  deep  and  shallow  water.  Sell  and  Hasselmann 
(1972)  have  shown  that  the  relative  values  of  the  fluxes  are  (+1),  (-4), 
and  (+3)  for  the  three  lobes  of  the  wave-wave  interactions  (Figure  5). 

Hence,  during  the  entire  time  that  the  wave  spectrum  is  adjusting  from 
-5  -3 

an  f  to  an  f  shape,  the  majority  of  the  momentum  is  going  into 
higher  frequencies  where  it  is  lost  to  turbulent  and  viscous  dissipa¬ 
tion.  If  the  deepwater  and  shallow-water  spectral  forms  are  integrated, 
the  total  energies  can  be  related  to  powers  of  the  peak  frequency.  The 
ratio  of  the  energies  of  the  spectra  can  be  defined  as 

r  =  !r  =  h  ™m2  (22) 

E  h  m 
s 

where  u>  is  the  radial  frequency  of  the  spectral  peak  (u)  =  2nf  ).  It 

m  ^  J  m  _^m 

should  be  understood  that  Equation  22  is  only  valid  when  the  f  as¬ 
ymptotic  spectral  shape  is  attained  in  shallow  water  |(h/g)^2U)m  <  l.oj; 
otherwise  the  more  general  form 

R  =  tanh2  k  h  (23) 

m 

*  Personal  communication,  C.  L.  Vincent  (1983),  U.  S.  Army  CERC,  Fort 
Belvoir,  Va. 
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Figure  5.  Schematic  representation  of  the  nonlinear 
wave-wave  interactions 

must  be  used.  Equations  22  and  23  both  assume  that  all  of  the  energy 
transferred  from  the  central  portion  of  the  spectrum  is  lost.  One- 
quarter  of  the  momentum  is  assumed  to  be  transferred  to  the  forward  face 
of  the  spectrum  where  it  should  not  be  immediately  lost.  Relating  the 
momentum  to  the  energy  in  the  spectrum  will  systematically  estimate  the 
amount  of  energy  lost  in  these  transformations.  In  shallow  water  the 
ratio  of  momentum  to  energy  becomes 


rE  OF  CHANGE  OF 
ENERGY  DENSITY 


where  M  is  momentum;  consequently  this  ratio  is  independent  of  fre¬ 
quency  and  one  should  expect  that  one-quarter  of  the  energy  transported 
out  of  the  central  portion  of  the  spectrum  is  added  to  the  energy  on  the 
forward  face.  Using  this  approximation,  the  net  energy  loss  due  to 
energy  fluxes  associated  with  wave-wave  interactions  is  given  by 


=  ![eo  - 


ET  =  7  E  -  tanh  (k  h)E 
L  4o  m  o 


and  the  remaining  energy  by 


■  Eo[i  ♦  l 


tanh  (kmh) 


34.  Equation  25  should  provide  a  good  approximation  to  the  energy 
loss  in  a  shoaling  spectrum  due  to  the  expected  change  in  spectral  shape 
as  long  as  the  spectral  peak  remains  at  about  the  same  frequency.  As 
the  spectral  peak  migrates  into  lower  frequencies,  the  region  of  energy 
loss  will  extend  into  these  frequencies  and  Equation  25  will  under¬ 
predict  the  amount  of  energy  loss.  Also,  it  should  be  borne  in  mind 
that  Equation  25  has  neglected  all  additional  energy  loss  mechanisms 
such  as  the  wave  bottom  interactions  mentioned  previously.  Thus  Equa¬ 
tion  25  is  expected  to  provide  a  reasonable,  slightly  conservative  ap¬ 
proximation  to  the  wave  energy  remaining  in  shallow  water.  Field  evi¬ 
dence  supporting  the  dominance  of  energy  loss  in  shallow  water  due  to 
this  single  mechanism,  along  with  a  more  detailed  theoretical  treatment, 
can  be  found  in  Vincent  (1981). 

35.  One  significant  unresolved  problem  still  remains  in  the  gen¬ 

eral  treatment  of  the  Phase  III  data,  due  primarily  to  the  assumption  of 
straight  parallel  offshore  contours.  As  illustrated  in  Figure  6, 
sheltering  effects  associated  with  the  geometry  of  the  coastline  can  in¬ 
troduce  significant  deviations  from  the  results  derived  from  a  straight 
coast  approximation.  The  solution  to  the  sheltering  problem  is  straight 
forward  for  both  sea  and  swell  wave  conditions.  Sea  wave  conditions  are 
described  by  a  two-dimensional  discrete  spectrum  (e^  A  function  can 
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Figure  6.  Generalization  of  spectral  and  wave  sheltering 

be  introduced  into  Equation  14  that  allows  only  certain  angle  bands  in 
the  two-dimensional  spectrum  to  arrive  into  a  sheltered  coastal  reach; 
thus  the  total  energy  can  be  expressed  as 

<26 
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where  A. .  is  defined  as 
ij 

a!  .  =  A. .A. 
iJ  iJ  J 

where  A.  =  1  for  wave  angles  that  can  arrive  at  the  coast  and  A.  =  0 
J  J 

for  wave  angles  that  are  prevented  from  reaching  the  sheltered  section 

of  the  coast.  Swell  wave  conditions  are  assumed  to  be  described  as  un¬ 
idirectional  monochromatic  waves;  thus  any  waves  propagating  from  a 
sheltered  area  are  set  to  zero. 

36.  The  amount  of  wave  sheltering  (or  the  selection  of  a  shelter 
ing  angle)  depends  on  the  location  of  the  input  Phase  II  station 
(Plates  1-6)  relative  to  the  Phase  III  station  being  considered.  Also, 
one  must  consider  the  location  of  the  land  boundary  points  within  the 
Phase  II  grid  relative  to  the  actual  shoreline.  There  is  a  limited 
amount  of  subjectivity  involved  in  the  selection  of  sheltering  angles 


when  all  of  the  above  variables  are  accounted  for. 


37.  In  very  shallow  water,  it  is  expected  that  the  similarity 

-3 

process  governing  the  equilibrium  range  k  form  must  eventually  break 
down  as  wave  breaking  begins  to  occur  at  frequencies  near  the  spectral 
peak.  At  that  point,  the  spectrum  may  evolve  to  a  different  shape  than 
that  given  by  Equation  21.  However,  evidence  from  a  study  by  Vincent 
(1981)  suggests  that  the  nearshore  energy  still  behaves  in  a  manner  con¬ 
sistent  with  the  equation.  The  integrated  form  of  Equation  21,  which 
predicted  a  dependence  of  the  wave  height  on  the  square  root  of  depth, 
provides  much  better  agreement  with  observed  wave  heights.  Therefore, 
in  shallow  water,  the  following  equation  can  be  used  to  describe  the 
variation  of  wave  height  with  depth 

H  =  f'1  (27) 

se  71  c 

max 


where  a  is  the  Phillip's  equilibrium  constant, 

f  =  0.9f  (28) 

c  m 

and  f^  is  defined  as  a  cutoff  frequency.  This  equation  may  be  used  to 

limit  the  locally  generated  sea  wave  height  found  in  the  Phase  III  wave 

data  without  having  to  rely  on  a  single  limiting  breaking  wave  condition 

This  has  the  advantage  that  all  of  the  wave  heights  exceeding  a  critical 

value  at  some  depth  are  not  mapped  into  the  same  nearshore  wave  height. 

1/2 

38.  Although  recent  field  evidence  appears  to  support  an  HMi 
relationship  in  very  shallow  water  for  broad  spectra,  there  is  a  sub¬ 
stantial  set  of  theoretical  and  laboratory  studies  that  support  a  linear 
relationship  between  monochromatic  wave  height  and  depth  (employed  for 
the  swell  breaking  limit), 


where  e  is  a  dimensionless  constant  with  a  value  generally  around  0.8. 
Thus  there  is  a  fundamental  difference  between  spectral  and  monochro¬ 
matic  forms  of  wave  breaking  in  shallow  water.  The  ability  to  treat  the 
local  sea  and  the  swell  wave  trains  separately  is  consequently  a  neces¬ 
sary  element  of  the  Phase  III  methodology. 

39.  Figure  7  presents  a  schematic  interpretation  of  the  overall 
Phase  III  methodology  as  developed  in  this  section.  Two  independently 
specified  wave  trains,  each  represented  by  an  appropriate  set  of  param¬ 
eters,  are  transformed  into  shallow  water.  This  separation  allows  the 
treatment  of  the  transformation  process  to  be  specified  differently  for 
each  wave  train.  Since  averaging  of  wave  approach  angles  is  not  always 
reasonable  and  since  different  mechanisms  appear  to  be  important  in 
broad  and  narrow  (monochromatic)  spectra,  this  approach  permits  the 
transformation  process  to  be  properly  posed  for  each  wave  train. 

Comparisons 

40.  Comparisons  to  actual  gage  measurements  are  necessary  to  de¬ 
termine  the  validity  of  the  long-term  wave  statistics  presented  in  WIS 
Report  9  (Jensen  1983).  Phase  III  wave  data  were  computed  for  166  sta¬ 
tion  locations  along  the  Atlantic  coast  (every  3  hr  for  the  20  years 
from  1956  through  1975).  Measured  wave  data  were  obtained  at  various 
locations  along  the  Atlantic  coast  from  the  wave  measurement  program  of 
the  Coastal  Engineering  Research  Center  (Thompson  1977). 

41.  Table  1  shows  the  latitude  and  longitude  of  the  locations  at 
which  measurements  are  available  and  the  corresponding  Phase  III  station 
locations.  Additional  information  concerning  the  type  of  wave  gage  used 
and  water  depth  at  the  gaging  site  is  included.  It  must  be  emphasized 
that  all  gage  measurements  are  highly  dependent  on  located  conditions 
such  as  water  depth,  bottom  topography,  temporal  changes  in  bottom  topog 
raphy,  and  pier  effects.  These  physical  conditions  are  considered  as 
subscale  features  in  comparison  with  the  Phase  III  shoreline  sections 
(10  miles*)  and  are  neglected  in  the  hindcast  data. 


To  convert  miles  to  kilometres,  multiply  by  1.609344 
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42.  The  Phase  III  wave  data  are  paired  as  closely  as  possible  to 
the  measured  data.  This  assumes  that  the  gage  data  (given  in  observa¬ 
tion  number  per  day  or  four  observations  per  day)  are  recorded  at  6-hr 
intervals  (0000,  0600,  1200,  1800  hours).  There  will  be  a  slight  phase 
shift  in  the  two  data  sets  caused  by  the  assumption  that  wave  conditions 
generated  via  the  Phase  III  methodology  neglect  propagation  time  from 
deep  to  shallow  water. 

WIS  Phase  III 
hindcast  wave  data 

43.  For  the  comparisons  in  this  report,  computed  wave  heights  are 

determined  from  the  combined  energy  associated  with  both  sea  and  swell 

[~~2  ^2 

wave  populations,  or  Hs  =*/H gEA  +  HSWELL  *  T^ie  Pea^  Peri°d  wave  infor¬ 
mation  presented  is  derived  from  the  selection  of  the  sea  or  swell  peak 
period  containing  the  dominant  wave  energy.  The  mean  directions  of  wave 
propagation  are  also  provided  for  sea  and  swell  wave  populations,  but  are 
not  used  because  no  measured  wave  data  contain  directional  information. 
Measured  wave  data 

44.  The  Coastal  Engineering  Research  Center  (and  its  predecessor, 
the  Beach  Erosion  Board)  has  operated  wave  gages  in  the  nearshore  area 
for  more  than  30  years.  Recent  wave  data  have  been  analyzed  and  tabu¬ 
lated  (Thompson  1977).  The  measured  gage  data  employed  in  all  com¬ 
parisons  are  derived  from  this  work.  Three  different  wave  recording  de¬ 
vices  (step-resistance,  continuous  wire,  and  pressure-sensitive)  were 
used  at  the  various  gage  sites  and  are  shown  in  Table  1.  The  approxi¬ 
mate  gage  locations  also  are  shown  in  Figure  8.  These  gages  have  been 
shown  to  respond  differently  to  identical  wave  conditions  (Thompson 
1977).  Step-resistance  wave  gage  results  were  found  to  be  approxi¬ 
mately  1  ft  higher  for  low  wave  conditions  and  20  percent  higher  for 
high  wave  conditions  in  comparison  with  continuous  wire  and  pressure- 
sensitive  gages.  This  will  produce  individual  bias  in  the  resulting 
wave  data . 

45.  The  wave  data  obtained  from  all  pressure-type  gages  are  based 
on  a  theoretically  derived  expression  relating  the  measured  dynamic  pres¬ 
sure  to  the  free  surface  oscillations.  The  basic  problem  encountered 


26 


was  how  to  compensate  for  the  pressure  signal  attenuation  due  to  gage 
submergence.  Generally,  this  relationship  is  based  on  linear  wave  theory 
where  a  simplified  transfer  function  acts  as  a  multiplier  to  change  a 
pressure  spectrum  to  a  free  surface  spectrum.  This  transformation 
caused  many  problems  in  the  analysis  of  the  wave  data  obtained  from  all 
pressure  gages  and  should  be  identified  at  the  onset  of  this  discussion. 
Due  to  its  sensitivity,  the  accuracy  of  the  Hg  derived  from  the  pres¬ 
sure  gages  may  be  questioned  (Thompson*). 

46.  All  summary  wave  data  (Hg  and  T  )  were  generated  from  digital 

wave  records  and  analyzed  in  a  similar  fashion.  The  significant  height 

is  equal  to  four  times  the  standard  deviation  of  record,  and  T  is 

p 

equal  to  the  frequency  at  ^lich  the  maximum  energy  density  exists.  The 
analysis  of  the  wave  record  employs  a  constant  frequency  interval  such 
that  converting  wave  frequencies  to  wave  periods  will  result  in  a  bias 
toward  8-  to  9-sec  peak  wave  periods  in  comparison  to  7-  to  8-  and  9- 
to  10-sec  peak  periods  (Thompson  1977). 

47.  Variations  in  water  depths  at  the  gage  location  also  will  in¬ 
fluence  the  measured  wave  conditions.  The  bottom  topography  in  a  near¬ 
shore  area  will  change  over  time  such  that  depth- limiting  wave  condi¬ 
tions  will  also  vary.  Changes  in  wave  conditions  from  10  m  (water  depth 
for  Phase  III  data)  to  7.6,  5.2,  4.6,  and  3.3  m  (water  depths  at  wave 
recorders,  as  shown  in  Table  1)  could  be  dramatic.  Wave  refraction  and 
shoaling  will  either  amplify  or  attenuate  wave  heights  as  much  as  25  per¬ 
cent.  Also,  the  nonlinear  energy  transfers  play  an  important  role  in 

the  nearshore  region,  affecting  both  the  H  and  T  result  (Vincent 

s  p 

1982).  One  also  must  consider  the  effects  of  wave  breaking,  on  a  wave- 
to-wave  basis,  within  a  given  wave  record  and  also  the  overall  collapse 
of  the  spectrum  in  very  shallow  water  or  when  extreme  wave  conditions 
prevail . 

48.  There  is  also  a  certain  degree  of  uncertainty  associated  with 
gage  records  containing  no  data  for  a  given  observation.  It  is  not 
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known  whether  a  malfunction  in  the  recorder  or  analysis  did  not  permit 

the  computation  of  the  wave  parameters  (H  and  T  ),  or  if  special  compu- 

s  p 

tations  were  not  made  because  the  significant  wave  heights  (derived  from 
evaluating  the  standard  deviation  of  the  wave  record)  were  less  than 
0.4  ft  (Thompson  1977).  This  uncertainty  would  lead  to  a  bias  in  the 
measured  data  toward  a  higher  percentage  of  large  wave  heights,  or 
omitting  wave  conditions  less  than  0.4  ft  (although  these  conditions 
were  considered  to  be  rare). 

49.  Thus  the  comparisons  presented  in  this  section  are  not  abso¬ 
lute  because  of  the  differences  in  physical  constraints  associated  with 
the  gaging  and  hindcasting  locations.  They  will  show  how  accurately 
(with  the  given  differences)  the  hindcast  data  represent  the  nearshore 
wave  climate. 

50.  Simple  yet  descriptive  techniques  were  employed  for  the  com¬ 
parisons:  time  series  plots,  cross  plots,  and  percent  occurrence  tables 
and  plots.  The  last  of  these  types  of  comparisons  is  probably  the  most 
indicative  of  the  ability  of  the  hindcast  wave  data  to  represent  a  wave 
climate  accurately.  The  two  data  sets  were  paired  according  to  the  gage 
observation.  These  paired-data  techniques  are  used  to  describe  the  ran¬ 
dom  error  of  the  hindcast  model.  However,  much  of  the  random  error  in 
these  computations  is  attributed  to  shifts  in  time  between  the  wave 
hindcasts  and  observed  wave  height.  Consequently,  for  the  evaluation  of 
these  hindcast  data  relative  to  wave  climate  representations,  major 
emphasis  is  centered  on  the  long-term  comparisons  of  means  and  probabil¬ 
ity  distributions  in  which  the  time  factor  has  been  removed. 

Time-history  comparisons 

51.  Although  the  time-history  comparisons  shown  in  Appendix  C  do 
not  provide  a  numerical  value  for  difference  or  similarity,  they  are 
qualitative  descriptors  of  the  random  error  of  the  hindcasting  technique. 
In  general,  all  of  the  time  series  plots  show  that  the  hindcast  Hg  is 
distributed  through  time  similar  to  the  distribution  of  the  measured 

Hg  (Figures  C1-C16) .  Some  intervals  show  the  computed  waves  to  be 
"out-of-phase"  with  respect  to  the  measured  data.  This  out-of-phase 
characteristic  can  be  attributed  to  the  Phase  III  methodology,  in  that 


propagation  of  wave  conditions  from  the  Phase  II  wave  data  to  the 
Phase  III  station  location  is  not  considered. 

52.  The  first  set  of  time-history  plots  (Figures  C1-C3)  represent 
comparisons  from  the  Nauset  gage  site.  Six  events  are  displayed  in 
these  figures  and  are  well  represented  by  the  hindcast  data.  Wave 
growth  over  time  in  the  hindcast  data  is  very  similar  to  that  displayed 
in  the  gage  results  for  all  six  events.  The  peak  wave  conditions  near 

9  December  1974  are  amplified  in  the  hindcast  results  in  comparison  with 
the  gage  data,  whereas  during  15-18  December  1974,  the  hindcast  data 
display  lower  maximum  conditions  over  the  gage  results.  The  reasons  for 
these  dissimilarities  can  be  attributed  to  the  gage  monitoring  highly 
site-specific  conditions  which  the  hindcasting  technique  was  not  de¬ 
signed  to  do. 

53.  The  second  set  of  time-history  plots  shows  comparisons  of 
gage  and  hindcast  results  at  Southhampton  (Figures  C4  and  C5).  There 
are  slight  phase  shifts  between  the  two  data  sets,  but  in  general  the 
hindcast  data  follow  the  wave  condition  trends  found  in  the  gage  re¬ 
sults.  Two  interesting  results  are  shown  in  these  plots.  The  first 
shows  that  the  gage  results  are  not  continuous.  Malfunctions  in  the 
recorder  or  analysis  cause  gaps  in  the  results.  For  the  five  storm 
events  shown  in  Figures  C4  and  C5,  only  one  storm  event  reaches  a  maxi¬ 
mum  condition  (25  November  1975).  When  data  are  missing  from  a  continu¬ 
ous  record  period,  a  certain  amount  of  bias  is  introduced  in  the  overall 
statistical  distributions  for  the  given  location.  The  second  point  that 
must  be  noted  is  that  during  periods  of  low  wave  energy  conditions,  the 
Phase  III  hindcast  data  will  become  very  calm  (i.e.  Hg  approaches  0.0), 
whereas  the  gage  records  larger  results.  These  differences  could  be  at¬ 
tributed  to  (a)  generation  of  wave  energy  in  the  nearshore  area  which 
the  gage  records  but  the  hindcast  assumes  to  be  negligible,  (b)  back¬ 
ground  noise  either  in  the  high  or  low  frequency  end  of  the  spectrum 
which  was  carried  through  the  analysis,  or  (c)  a  result  of  the  physical 
geometry  of  the  given  gage  site. 

54.  The  third  set  of  time-history  plots  presents  comparisons  of 
gage  and  hindcast  results  at  Nags  Head  (Figures  C6-C8) .  The  peak  wave 


periods  (T  =  1/f  )  derived  from  the  gage  and  hindcast  results  also  ai?e 
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displayed.  The  results  clearly  show  that  the  hindcast  accurately 

represents  storm  event  conditions.  The  growth  and  decay  of  the  near¬ 
shore  wave  conditions  are  also  adequately  described.  The  peak  period 
results  display  dissimilar  trends  between  the  gage  and  hindcast  data 
sets.  The  primary  reason  for  the  differences  are  caused  by  the  gage 
recording  a  multicomponent  (in  frequency)  wave  system,  whereas  the  hind¬ 
cast  data  are  represented  by  a  two-component  system  (sea  and  swell). 
Also,  the  one-dimensional  frequency  spectrum  derived  from  the  gage 
record  contains  89  frequency  bands  (at  a  constant  interval),  whereas 
the  hindcast  spectrum  is  modeled  with  20  frequency  bands  (at  variable 
frequency  intervals)  (Resio  1981).  Thus  the  declaration  of  the  peak 
period  results  obtained  from  the  gage  analysis  would  show  larger  changes 
from  observation  to  observation. 

55.  Three  additional  sets  of  time-history  Hg  results  are  shown 
in  Figures  C9-C16  for  Wrightsville  Beach,  Holden  Beach,  and  Daytona 
Beach.  The  trends  in  the  hindcast  data  are  in  agreement  with  the  gage 
results.  Because  of  the  difference  in  water  depths  (10  m  for  the  hind¬ 
cast  data  and  5.18,  4.57,  and  3.35  m  for  Wrightsville,  Holden,  and  Day¬ 
tona  Beach  data,  respectively),  the  hindcast  H  for  maximum  conditions 

s 

generally  overpredicts  the  gage  data  because  of  depth-limited  wave  break 
ing,  whereas  during  moderate  to  low  wave  conditions  the  differences  be¬ 
tween  the  two  data  sets  are  minimal.  Even  slight  variations  in  the  wave 
climate  represented  in  the  gage  record  (Figure  C14,  from  Daytona  Beach) 
are  very  closely  approximated  in  hindcast  data. 

Percent  occurrence  comparisons 

56.  There  is  a  degree  of  uncertainty  associated  with  the  statis¬ 
tical  representations  of  the  wave  climate  at  a  given  location.  One  must 
be  aware  that  even  the  gage  cannot  uniquely  describe  actual  wave  condi¬ 
tions  during  a  given  sampling  period.  Use  of  long-term,  continuous  (de¬ 
fined  here  as  a  constant  number  of  observations  per  day  over  a  prolonged 


period  of  time)  Hg  and  T^  data  will  offset  the  inconsistency  in  the 
actual  and  measured  wave  conditions.  Unfortunately,  the  longest  mea¬ 
sured  wave  data  set  is  from  Nags  Head,  North  Carolina  (spanning  7  years 


with  a  total  of  4,492  H  and  T  observations).  This  means  that 
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4,268  observations  are  either  calm  (defined  by  Thompson  1977)  or  missing 
The  working  data  set  only  contains  approximately  one-half  of  the  total 
expected  number  of  observations.  Thus  a  bias  of  unknown  quantity  is  in¬ 
troduced  into  the  gage  results  because  of  the  noncontinuous  data  set. 

The  hindcast  data  presented  in  the  percent  occurrence  plots  also  will  be 
biased  because  they  are  paired  as  closely  as  possible  to  the  gage  data 
through  time.  Hence,  any  differences  found  in  these  comparisons  will  be 
caused  by  the  hindcast' s  own  bias.  The  percent  occurrence  in  the  Hg 
data  will  also  be  influenced  by  the  water  depth  differences  between  the 
hindcast  results  (10  m)  and  the  gage  results  (9.76,  5.18,  4.57,  and 
3.35  m).  All  wave  conditions  will  be  subject  to  the  effects  of  wave 
shoaling  from  10  m  to  the  gage  water  depth  where  increases  as  much  as 
25  percent  can  be  expected  based  on  theoretical  considerations.  Wave 
refraction  can  diminish  the  effects  of  wave  shoaling,  but  under  certain 
conditions  an  increase  in  wave  height  can  be  expected. 

57.  For  each  site,  the  percent  occurrence  for  and  T^  is 

compared.  The  paired  data  (hindcast  and  gage  results)  are  separated 

into  H  intervals  of  0.5  m  and  the  T  in  intervals  of  1.0  sec. 
s  p 

These  plots  are  generated  for  comparison  purposes  only  and  should  not  be 
construed  as  the  representative  wave  climate  expected  for  that  specified 
location. 

58.  Figure  9  displays  the  percent  occurrence  of  Hg  and  T^ 
conditions  for  Nauset  Beach,  Massachusetts,  from  1974-1975.  The  wave 
recorder  is  a  pressure  gage  mounted  approximately  10  m  below  the  mean 
water  level.  The  hindcast  Hg  results  are  shown  to  underpredict  the 
gage  results  in  all  but  the  first  category.  The  difference  in  the  per¬ 
cent  occurrence  for  large  wave  conditions  could  be  attributed  to  the 
greater  water  depth  found  at  the  gage  location.  The  reason  for  the 
large  difference  in  the  first  two  categories  (<0.5  m  and  <1.0  m)  cannot 
be  explained  at  this  time.  The  tail  of  the  percent  occurrence  (for  wave 
conditions  from  <2.0  to  <5.0  m)  shows  that  the  hindcast  varies  from  the 
gage  results  by  4  percent  at  its  maximum.  The  percent  occurrence  in  the 
measured  T  results  demonstrates  that  there  is  a  bias  away  from  wave 


periods  in  the  category  <8  sec.  There  is  also  a  large  portion  (45  per¬ 
cent)  of  wave  conditions  with  T  values  <4  sec.  By  linear  theory, 
wave  frequencies  greater  than  0.27  Hz  should  not  be  resolvable  in  a 
water  depth  of  10  m.  There  appears  to  be  some  question  as  to  the 

validity  of  all  results  with  T  values  less  than  3.7  sec.  There  were 

P 

problems  associated  with  initial  measurements  taken  at  Nauset  Beach  (as 
previously  mentioned,  paragraph  45)  that  would  result  in  inconsistent 
wave  conditions.  Rather  than  edit  the  measured  data  by  theoretically 
derived  limiting  conditions,  it  was  felt  that  the  gage  results  should  be 
presented  in  untampered  form. 

59.  The  second  location  for  comparisons  is  the  Southhampton,  New 
York,  gage  site.  The  gage  is  a  pressure  type  gage  located  in  approxi¬ 
mately  10  m  of  water.  The  percent  occurrence  of  Hg  and  T^  is  shown 
in  Figure  10.  The  hindcast  results  compare  extremely  well  with  the  gage 
results  in  terms  of  the  Hs  data.  The  maximum  difference  is  approxi¬ 
mately  8  percent  while,  in  general,  differences  of  approximately  2  per¬ 
cent  exist  throughout  the  distribution.  The  maximum  difference  in  the 
Tp  results  is  approximately  10  percent  found  in  the  gage  biased 
category  <8  sec. 

60.  The  third  gage  location  for  comparisons  is  at  Nags  Head, 

North  Carolina  (Figure  11).  Two  types  of  monitoring  devices  were  em¬ 
ployed.  A  step-resistance  gage  (from  1968-1972)  and  a  continuous  wire 
gage  (from  1972-1975).  The  differences  found  in  the  percent  occurrence 
plot  for  the  Hg  categories  <0.5  and  <1.0  m  can  be  attributed  to  the 
bias  associated  with  the  step-resistance  gage  (measuring  approximately 
0.3  m  high  for  low  wave  conditions).  If  this  "error"  were  used  to  com¬ 
pensate  the  gage  data,  then  the  differences  between  the  hindcast  and 
gage  data  would  be  approximately  10  percent.  The  percent  occurrence 
distribution  for  larger  wave  conditions  shows  differences  of  5  percent 
(at  the  maximum)  between  the  two  data  sets.  Because  of  the  differences 
in  the  water  depth  where  the  wave  information  is  generated  (10  m  for  the 
hindcast  results  and  5.18  m  for  the  gage  results),  there  will  be  slight 
dissimilarities  in  the  extremely  large  wave  heights.  One  possible  cause 
for  the  difference  could  be  the  difference  in  the  techniques  used  in  the 
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Figure  11.  Percent  occurrence  versus  H  and  T  for 
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time-paired  (4492  observations)  measured  and  computed 
data  sets  at  Nag  Head,  North  Carolina 


assignment  of  peak  period  in  the  gage  hindcast  analyses.  In  Figure  C6, 
the  gage  records  indicate  a  shift  in  period  from  ~3  sec  to  M4  soe  while 
the  wave  height  remained  almost  constant.  Shifts  in  wave  per. ,d,  as 
those  shown  in  Figure  C6,  are  considered  to  be  unrealistic,  and  for  wave 
heights  of  0.2  to  0.4  m,  relatively  unimportant.  In  Figure  C8,  t^e  ..'age 
recorded  period  shifts  from  VJ  sec  to  "V4  sec  within  one  recorded  inter¬ 
val  while  the  wave  height  was  increasing  slightly.  Again,  this  shift  is 


31 


not  realistic;  and  the  hindcast  representation  is  more  appropriate. 

Also,  the  differences  in  the  results  could  be  caused  by  the  input, 

deepwater  wave  information.  The  hindcast  model  (Resio  1981)  tends  to 
"wash-out"  long-period  low-amplitude  swell  where  the  energy  is  numeri¬ 
cally  dispersed  throughout  the  grid. 

61.  The  statistical  distribution  comparisons  for  H  and  T  at 

s  p 

Wrightsville  Beach,  North  Carolina,  are  presented  in  Figure  12.  At  this 

location,  two  wave  gaging  systems  were  employed  (a  step-resistance  type, 

from  1970  to  1974,  and  a  continuous  wire  type,  from  1975  thereon). 

Again  the  two  data  sets  compare  favorably  for  H  conditions  from 

s 

<1.5  m  to  their  extreme,  while  the  bias  in  the  step-resistance  gage 
(1.0  ft  for  low  wave  conditions)  produces  the  differences  found  in  the 
<0.5  and  <1.0  m  height  increments.  If  the  bias  in  the  gage  data  is  cor¬ 
rected,  the  difference  between  the  two  data  sets  is  20  percent  in  wave- 
height  classes  <0.5  and  <1.0  m.  Although  this  difference  is  very  large, 
the  total  percent  occurrence  of  wave  conditions  in  both  categories  dif¬ 
fers  by  only  4  percent.  Similar  inconsistencies  are  found  in  the  T 

P 

distributions  as  were  previously  found  in  the  Nags  Head  data. 

62.  The  Holden  Beach,  North  Carolina,  and  Daytona  Beach,  Florida, 

percent  occurrence  plots  are  displayed  in  Figures  13  and  14.  The  water 

depths  at  the  gage  location  are  4.57  and  3.35  m,  respectively,  and  the 

wave  monitoring  devices  for  these  locations  were  continuous  wire  type 

(Holden  Beach)  and  pressure  type  (Daytona  Beach).  Again,  under  low  wave 

conditions  the  hindcast  produces  a  larger  percentage  of  Hg  <0.5  m 

while  the  gage  results  show  a  larger  percentage  of  Hg  <1.0  m.  The 

differences  between  the  two  sets  of  data  are  approximately  15  percent  at 

their  maximum.  Comparing  the  percent  occurrence  of  T  derived  from 

P 

the  hindcast  and  gage  data  one  finds  good  agreement  at  Holden  Beach  for 
nearly  all  wave  period  categories  (with  exception  of  the  biased  <8  sec 
category  of  the  gage  data).  At  Daytona  Beach  the  gage  results  show  a 
higher  percentage  of  long  wave  periods  whereas  the  hindcast  data  consist 
of  a  higher  percentage  of  shorter  period  wave  conditions. 

Extreme  storm  event  comparisons 


63.  A  final  comparison  is  made  between  extreme  storm  wave 


WRIGHTSVILLE  BEACH  GAGE  LOCATION 
1971-1975 

LEGEND 

- GAGE 

— - HINDCAST 


GAGE  AND  HINDCAST  ZERO  PERCENT 


Tp,  SEC 


Figure  12.  Percent  occurrence  versus  H  and  T  f 
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time-paired  (4126  observations)  measured  and  computed 
data  sets  at  Wrightsville  Beach,  North  Carolina 


conditions.  The  data  were  derived  by  selecting  the  top  20  to  30  storm 

events  within  the  gage  sampling  period.  The  maximum  Hg  from  that 

event  was  recovered .  The  date  of  these  events  then  was  used  to  scan  the 

hindcast  data  set  to  select  the  maximum  H  from  the  storm.  These  max- 

s 

imum  Hg  conditions  are  not  identically  paired  in  time  (lag  time  caused 

by  nonpropagation  of  wave  conditions  from  the  Phase  II  input  conditions) 

but  are  representative  of  the  event.  In  some  instances,  the  water  depth 

at  the  gage  site  will  greatly  influence  the  maximum  wave  conditions. 

These  data  will  be  biased  toward  lower  H  results  for  a  particular 

s 

storm  condition  in  comparison  to  the  hindcast  results. 

64.  The  Nauset,  Massachusetts,  gage  site  extreme  Hg  conditions 
are  compared  with  the  hindcast  results  shown  in  Figure  15.  For 
conditions  up  to  approximately  2.5  m,  the  hindcast  results  generally  tend 
to  be  slightly  underpredictive  to  the  gage  data.  Above  this  point,  the 
gage  results  are  shown  to  be  from  1.0  to  2.0  m  higher  for  maximum  wave 


Figure  15.  Cross  plot  of  measured  (h  =  10.7  m)  and  computed 
(h  *  10  m)  maximum  H  conditions  occurring  in  a  specific 

event  for  Nauset,  Mass. 


conditions  in  comparison  with  the  hindcast  results.  From  this  data  set, 
only  one  storm  event  was  completely  missed  by  the  hindcast  (H  sr  2.0 

and  H  2  0.6  1);  all  other  events  were  found  in  the  h^n§cast 

.  *  SHINDCAST 
data  set . 

65.  The  plot  of  extreme  storm  conditions  for  Southhampton  is  dis¬ 
played  in  Figure  16.  These  data  show  that  the  hindcast  results  are 
underpredictive  for  low  wave  conditions  (less  than  approximately  2.0  m) , 

then  tending  toward  over predicting  H  results  greater  than  2.5  m. 

s 

There  is  no  apparent  reason  for  these  differences. 

66.  The  final  four  figures  (Figures  17-20)  are  comparisons  of  the 
extreme  wave  conditions  for  Nags  Head,  Wrightsville  Beach,  Holden  Beach, 
and  Daytona  Beach.  All  gage  results  are  highly  dependent  on  the  water 
depth  at  the  gage  sites  (5.18,  5.18,  4.57,  and  3.35  m,  respectively) 

and  become  biased  toward  lower  H^  values  in  comparison  with  the  hind¬ 
cast  results  and  tend  to  underpredict  the  gage  results.  Above  this 
point,  the  hindcast  results  (being  in  deeper  water)  tend  to  overpredict 
the  measured  data.  In  these  comparisons  there  is  only  a  limited  number 
of  storm  events  completely  missed  by  the  hindcast,  the  maximum  number 
occurring  at  the  Wrightsville  Beach  gage  site. 

Summary  and  Conclusions 


67.  The  methodology  of  the  Phase  III  portion  of  the  Wave  Informa 
tion  Study  has  been  summarized  and  comparisons  have  been  made  to  six 
unique  gage  locations  along  the  Atlantic  coast.  The  transformation  ma¬ 
trix  operator  generating  the  hindcast  wave  data  remained  constant  from 
one  gaging  site  to  the  next.  Only  variations  in  the  first-order  approx 
imation  to  the  shoreline  orientation  (wave  generating  model)  slightly 
change  the  generalized  transformation  matrix  operator. 

68.  There  are  areas  along  the  Atlantic  coast  where  the  Phase  III 
wave  data  (WIS  Report  9,  Jensen  1983)  may  not  be  appropriate.  Since 
wave  conditions  are  not  measured  along  the  total  length  of  the  Atlantic 
coast,  all  Phase  III  station  data  cannot  be  verified.  The  whole  coast 
of  Maine  presents  a  significant  problem  in  the  declaration  of  10-mile 
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Figure  20.  Cross  plot  of  measured  (h  =  3.4  m)  and  computed 
(h  =  10  m)  maximum  H  conditions  occurring  in  a  specific 
event  for  Daytona  Beach,  Florida 

straight  shoreline  segments.  This  coastline  has  numerous  bays  and  off¬ 
shore  islands  which  negate  the  straight  and  parallel  bottom  contours 


assumption.  For  this  area,  all  Phase  ITT  shoreline  station  segments  are 
aligned  with  the  outermost  offshore  islands.  Embayments,  tidal  inlets, 
and  sounds  are  neglected  whenever  they  exist  within  a  given  Phase  III 
station  section.  In  some  instances,  adjustments  must  be  made  for  large 
offshore  shoals  (for  example,  Nantucket  Shoal,  Massachusetts,  and  Frying 
Pan  Shoal,  North  Carolina).  Also,  wave  conditions  west  of  Cape  Cod, 
Massachusetts,  may  only  be  a  gross  approximation  of  the  actual  wave  con¬ 
ditions.  Table  2  lists  certain  Phase  III  station  locations  (shown  in 
Plates  1-6)  along  the  Atlantic  coast  where  "site-specific"  modeling 
should  be  utilized  to  accompany  the  Phase  III  wave  information. 

69.  The  comparisons  between  gage  and  hindcast  results  demon¬ 
strated  that  the  accuracy  of  the  Phase  III  hindcasting  technique  is  not 
lost  while  increasing  computational  efficiency.  These  comparisons  have 
shown  that  the  hindcast  represents  existing  wave  conditions  (time-paired 
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and  statistically)  considering  all  variables  associated  with  gaging 
operations.  As  in  the  case  of  the  Phase  I  hindcast  wave  data  (Corson 
and  Resio  1981),  the  shallow-water  (Phase  III)  wave  data  did  not  consis 
tently  overpredict  nor  underpredict  the  measured  Hg  conditions,  al¬ 
though  the  hindcast  data  showed  a  slight  bias  toward  lower  wave 

periods.  The  absolute  cause  of  these  differences  cannot  be  defined; 
therefore  a  clear-cut  adjustment  to  the  hindcast  data  cannot  be 

made.  The  hindcast  wave  data  presented  in  WIS  Report  9  (Jensen  1983) 
should  be  used  "as  is"  for  ocean  engineering  projects. 
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Table  1 

Physical  Parameters  Associated  with  Gage  Sites 


Gage  Site 
Name 

Location 

Data  Set 
Dates 

Water 

Depth 

m 

Gage  Type 

Nauset,  Mass. 

41.8°N,  69 . 93°W 
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10.67 

Pressure 

(1/74-12/75) 

Southhampton , 
N.  Y. 

40.85°N,  72.42°W 

1973-1975 

9.76 

Pressure 

(N/A) 

Nags  Head 

N.  C. 

35.92°N,  75.60°W 

1968-1975 

5.18 

Step- resistance 
(19/63-10/72) 
Continuous  wire 
(11/72-12/75) 

Wrightsville 
Beach,  N.  C. 

34. 22°N,  77 . 78°W 

1971-1975 

5.18 

Step- resistance 
(3/70-12/74) 
Continuous  wire 
(1/75-12/75) 

Holden  Beach, 
N.  C. 

33.92°N,  78.30°W 

1971-1975 

4.57 

Continuous  wire 
(2/71-12/75) 

Daytona  Beach 
Fla. 

29 . 15°N,  80 . 97°W 
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Plates  1-6  for  geographical  location  of  the  Phase  III  stations. 


Locations  of  Phase  III  stations  for  shallow-water  wave  information 
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Plate  2.  Locations  of  Phase  III  stations  for  shallow-water  wave  information 

along  the  Atlantic  coast,  Region  2 
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Plate  3.  Locations  of  Phase  III  stations  for  shallow-wave  wave 
information  along  the  Atlantic  coast,  Region  3 
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Appendix  A:  Mathematical  Treatment  of  the  Transformation  of 
Probabilities  from  One  Coordinate  System  into  Another 


1.  Given  an  operator  such  that  one  set  of  parameters  ,  ...  , 
can  be  uniquely  mapped  into  another  coordinate  system  ,  ...  ,  , 


W!  =  ,  ...  ,  W  ) 

1  Ti  1*  n 


wi  =  *  Wn} 


then  a  probability  density  function  in  W  can  be  transformed  into  a 
probability  density  function  in  W'  , 


p(w*,  ...  ,  w;) 


=  p[w  •••  •  V . VW1 . v]  |Ji 


where  J  is  the  Jacobian  of  the  transformation 
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For  an  increment  of  probability,  the  relationship  between  probabilities 
in  the  different  parameter  spaces  can  be  written  as 


•  •*>! 


vV 


Or,  in  terms  of  mapping  a  continuous  function  from  W  into  W'  would 
yield 


V  ) 
n 


dw: 


dW 

n 


hi 


p(w 


i’ 


V 


dW, 


dW 


(A5) 


where  the  limits  of  region  R  are  defined  by  applying  the  operator  <|> 
to  the  limits  of  S  .  Vice  versa,  the  limits  of  S  can  be  obtained  by 
applying  the  inverse  operator  4*  ^  to  R  • 

2.  Numerical  estimation  techniques  for  probabilities  generally 
deal  with  a  discretized  representation  of  the  probability  density  func¬ 
tion;  consequently,  one  needs  to  transform  Equation  5  in  terms  of  a  dis¬ 
cretized  approximation.  Given  r  as  an  incremental  index  to  the  loca¬ 
tion  of  a  discrete  category  in  W,  probability  increments  of  r  can 
then  be  mapped  into  probability  increments  in  r'  in  an  analogous 
manner  to  Equation  5.  Letting  p'  denote  the  discretized  probability 
density  function  in  r'  space, 


pUj, 


rn)  (A6) 


where  (j^,  k^),  ...  ,  ( j n ,  k^)  are  the  discretized  limits  in  region  S 
in  W  corresponding  to  a  given  incremented  region  R  in  W’  . 

3.  Thus  when  a  small  set  of  parameters  is  sufficient  to  cate¬ 
gorize  a  probability  function  for  practical  purposes,  it  is  feasible  to 
approach  the  transformation  of  probabilities  from  one  such  parameter 


space  into  a  related  parameter  space  by  storing  indices  in  one  coordi¬ 
nate  system  as  a  function  of  indices  in  the  related  coordinate  system. 
This  method  is  particularly  advantageous  when  the  form  of  <J>  is  not  a 
simple  analytic  function  of  all  parameters.  In  such  a  case,  the  evalua¬ 
tion  of  the  Jacobian  must  be  done  numerically,  which  can  be  exceedingly 
tedious.  Such  complex,  nonlinear  forms  for  <{>  are  typical  in  nearshore 
wave  transformation  processes.  An  approximation  of  the  following  form 
can  be  used  for  such  transformations. 

4.  First,  determine  rj,  =  (r^,  ...  r  )  This  represents  a  map¬ 
ping  of  the  central  point  of  the  discretized  category  r^,  • • •  rn  into 
a  discretized  category  rj^  ...  r^  .  For  each  discretized  category  in 
parameter  space  r  there  will  be  a  unique  corresponding  discretized 
category  in  r'  .  In  this  manner,  an  approximation  for  a  probability 
density  function  in  r '  can  be  formed  from  a  probability  density  func¬ 
tion  in  r 

p(rj,  •  •  •  ’  rn^  =  ^2  '  '  '  P^rl’  *  ’  *  ’  rn^6  l*K£)  ‘  r'  ]  (A7) 

r,  r 

1  n 

where  the  delta  function  6[4>(r)  -  r* ]  ensures  that  ...  dp  is  con¬ 
served  in  the  mapping.  Since  Equation  7  represents  the  mapping  of  dis¬ 
crete  probability  increments  rather  than  probability  densities,  some 
bias  will  be  generated  by  this  operation;  however,  if  the  incrementiza- 
tion  is  small  the  bias  should  be  minimal.  Also,  additional  accuracy  can 
be  obtained  by  modifications  to  Equation  7  which  consider  a  gradient  in 
p  across  an  increment  in  r'  .  For  the  present,  Equation  7  will  be 
used  as  written  to  obtain  the  approximation  to  p(r')  ,  with  the  excep¬ 
tion  that  in  practice  the  delta  function  behavior  is  obtained  numeri¬ 
cally  by  stepping  through  indices  r  (only  once)  into  its  appropriate 
category  in  r'  . 


Appendix  B:  Comparison  of  Spectral  and  Monochromatic 


Representation  for  Refraction  and  Shoalin 


1.  Since  most  engineers  are  more  familiar  with  the  monochromatic 
representations  for  refraction  and  shoaling,  a  comparison  with  the  spec¬ 
tral  representations  is  included  here,  along  with  a  comparison  of  the 
differences  predicted  by  the  two  methods.  The  traditional  conservation 
of  energy  flux  in  monochromatic  wave  is  written  as 


(cgEb)  = 


[c  Eb] 

V8  b 


where  c  is  the  group  velocity,  E  is  wave  energy,  b  is  the  width 
between  two  adjacent  wave  rays,  and  the  subscript  o  refers  to  initial 
conditions  (deep  water).  Figure  B1  shows  two  wave  rays  coming  into  a 
coast.  The  spacing  between  the  rays  can  be  shown  to  be  equal  to  a  pro¬ 
jection  of  an  orthogonal  onto  a  plane  parallel  to  the  depth  contours; 


hence  the  ratio  of  b 
the  angles  6  and  6 


can  be  defined  in  terms  of  the  cosines  of 


cos  6 
cos  6 


where  6q  and  6  are  the  initial  angle  and  subsequent  angles  along  the 
wave  ray,  respectively.  This  ratio  bQ/b  is  commonly  referred  to  as 
the  refraction  coefficient  and  the  ratio  c  /c  ,  the  shoaling  coef- 

8q  8 

ficient.  From  Equation  B1  and  B2,  the  energy  along  the  ray  can  be 


written  as 


2.  If  we  now  examine  a  narrow  spectrum  approaching  in  the  limit  a 
delta  function  in  angle  and  frequency.  Equation  12  yields 
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(B4) 


where  E2(f,6)  is  the  initial  energy  distribution  given  by 

E2(f,0)  =  Eo6(0-0’)6(f-f’) 


where  0'  is  the  central  angle  of  the  angular  distribution  and  f'  is 
the  frequency  of  the  spectral  peak.  Integrating  over  both  delta  func¬ 
tions  gives 


E  =  E 
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36 


36 


(B5) 


From  Snell's  law  we  have: 


sm  0  =  —  sm  0 
c  o 

o 


(B6) 


Differentiating  both  sides  with  respect  to  0  gives 


cos  0  =  — -  cos  0 
c 

o 


30 


30 


(B7) 


which  when  rearranged  provides  an  evaluation  for  the  partial  derivative 
in  Equation  B5. 
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(B8) 


Substituting  this  into  Equation  B5  gives 


E  =  E 


cos  0  g 
o  o 


o  cos  0  c 


(B9) 


g 


which  is  identical  with  Equation  B3,  as  it  should  be  for  this  case. 
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Figure  C7.  Plots  of  Hg  and  versus  time  for  a  period  in 

December  1971  for  Nags  Head,  North  Carolina 
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versus  time  for  a  period  in  September-October  1974  for 
Holden  Beach,  North  Carolina 
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Appendix  D:  Notation 


Phase  velocity  o£  the  wave  component 
Group  velocity  of  the  wave  component 
Total  deepwater  spectral  energy 
Total  shallow  water  spectral  energy 
Total  integrated  spectral  wave  energy 

Total  spectral  energy  loss  due  to  wave-wave  interactions 

One-dimensional  spectrum 

Two-dimensional  spectrum 

Frequency 

Cutoff  frequency 

Frequency  of  the  spectral  peak 

Two-dimensional  wave  spectrum  in  wave  number  (k)  space 

Acceleration  due  to  gravity 

Water  depth 

Wave  height 

Frequency  increment 

Incremental  values  of  the  deepwater  parameter 
Angle  increment 
Wave  number  modulus 
Wave  number  space 

Number  of  intervals  used  to  resolve  the  i*"*1  parameter 
Momentum 

Number  of  parameters 

Action  density  defined  in  terms  of  the  wave  number 
modulus  k 

Total  number  of  transformations 
Matrix  form  of  the  transformation  processes 
Nearshore  wave  parameter  value 
Nondimens ional  dispersion  relationship 
Wave  period 

Transformation  matrices 
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Definition 


Symbol 
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e 

0 

o 

6 

X 

A.  . 
ij 

4> 

U) 

0) 

Subscript 


Peak  spectral  wave  period 
Phillips  equilibrium  constant 
Spectral  shape  function 
Coefficient  flagging  wave  sheltering 

Dimensionless  constant,  dependent  on  the  breaking  limit 
Direction 

Deepwater  wave  angle 

Central  angle  of  the  spectrum 

Constant 

Matrix  of  precomputed  values  of  p  (f,h)  in  discrete  form 
Normalized  angular  distribution  function 
Nondimensional  function 
Radial  frequency  of  spectral  peak 


o 


Deepwater  conditions 


